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Abstract An exact solution of Einstein's equations representing the static 
gravitational field of a quasi-spherical source endowed with both mass and 
mass quadrupole moment is considered. It belongs to the Weyl class of so- 
^ ' lutions and reduces to the Schwarzschild solution when the quadrupole mo- 

ment vanishes. The geometric properties of timelike circular orbits (includ- 
ing geodesies) in this spacetime are investigated. Moreover, a comparison 
between geodesic motion in the spacetime of a quasi-spherical source and 

■ non-geodesic motion of an extended body also endowed with both mass and 
mass quadrupole moment as described by Dixon's model in the gravitational 
field of a Schwarzschild black hole is discussed. Certain "reciprocity relations" 
between the source and the particle parameters are obtained, providing a 
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further argument in favor of the acceptability of Dixon's model for extended 
bodies in general relativity. 

Keywords Extended bodies in general relativity ■ Dixon's model 
PACS 04.20.Cv 



1 Introduction 

The fully relativistic multipole moments of a stationary spacetime have been 
introduced by Hansen [I], generalizing previous results by Geroch [5] valid 
for static spacetimes only. Hansen's formulation reduces to Geroch's one in 
the static limit, in the sense that the recursive definitions of moments are the 
same but with a different potential. Beig [3] contributed to clarify Hansen's 
approach considering a different definition of center of mass; in this way 
the expansion of the Hansen moments around the center of mass deter- 
mines the multipole moments uniquely. Beig and Simon [4j[5] also applied 
the above mentioned definition to the case of stationary axisymmetric space- 
times. Thorne [6J, before the works of Beig and collaborators, gave another 
definition of multipole moments which later Gursel [7] has shown to be equiv- 
alent to Hansen's definition for a source with nonzero rest mass. Many exact 
solutions of Einstein's equations for sources having multipolar structure as 
described by the Geroch-Hansen formulation are known, mostly belonging to 
the Weyl class of stationary axisymmetric spacetimes 8 and obtained with 
a suitable use of Ernst potentials and generating techniques [9], and hence 
formally very complicated. 

On the other hand, the motion of extended bodies (considered as test 
bodies, i.e. with backreaction neglected) in any given background was well 
established after the works of Mathisson, Papapetrou (up to the dipolar 
structure) pIJIfTT] and Dixon (including any multipolar structure) 12,13,14, 

EEMU]- 

In this paper we analyze circular orbits (including geodesies) in the grav- 
itational field of a body endowed with both mass and mass quadrupole 
moment, which reduces to the familiar Schwarzschild solution when the 
quadrupole moment vanishes. Then we compare geodesic motion of a test 
particle on the equatorial plane with that of an extended body also endowed 
with both mass and mass quadrupole moment as described by Dixon's model 
in the gravitational field of a Schwarzschild black hole. We investigate the 
correspondence between the source and the particle parameters. We obtain 
certain "reciprocity relations" leading to the identification of Dixon's model 
quadrupole parameters with those underlying Geroch-Hansen approach. This 
is a novel result which should be regarded as a further argument in favour of 
the acceptability of Dixon's model for extended bodies in general relativity. 

2 Test particles in the field of a quasi-spherical source 

The metric of a nonrotating mass with a quadrupole moment has been ob- 
tained long ago by Erez and Rosen [17j , later corrected for several numerical 
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coefficients independently by Doroshkevich et al. [18] and Young and Coul- 
ter [19] . It is a solution of the static Weyl class of solution with the metric 
element of the following form 

ds 2 = -e^dt 2 + e 2 ^\dp 2 + dz 2 ) + p 2 e~ 2 ^dz 2 , (1) 

with i\) and 7 functions of p and z only. The associated vacuum Einstein's 
equations are 

Vw + -tp P + 4>zz = , 

7p = p(4> 2 p -4>l), j z = 2pip p ip z . (2) 

Within this class, the solution representing a source with mass (monopole) 
and quadrupole structure can be written in the following form by using pro- 
late spheroidal coordinates x and y 



^ = <Ao + q\($y 2 - 1) 



i(3x 2 -l)Vo + ^ 



(3) 



7 = 70 + q [2 7o - 3(1 - y 2 ) (xfa + 1)] 

+q 2 1 To + ^(1 - y 2 ) [3(x 2 - l) 2 ^ 2 + 2x(3x 2 - 5)^o + (3x 2 - 4)] 

~ ~ y2) [ (1 " 3:2)2(1 _ 9x2)2 ^ 2 ° + 2x{9x2 " 7)v, ° + {3x2 ~ 4) ] 

where 

are the monopole quantities corresponding to the Schwarzschild solution 
(which is obtained for q = 0). The relation with Weyl coordinates p and 
z is given by 

x = (r_L + r_) , w = (r + — r_) , (5) 

with 

r ± = [p 2 + (z±M) 2 ] 1 / 2 . (6) 

Note that on the symmetry hyperplane z — we have r + = r_ = [p 2 + 
M 2 ] 1 / 2 = r* and hence x = r*/M, y = 0. 

It is worth to mention that there exist in the literature several different 
solutions corresponding to a quasi-spherical source, all with the correct grav- 
itational potential of a massive static source with a quadrupole moment in 
the Newtonian limit [201121]. For these metrics some properties of geodesic 
motion have been analyzed [22ll23| . Geodesies in the Erez-Rosen spacetime 
have been also extensively investigated [2"4"1[2"5] . A generalization to the sta- 
tionary case including rotation of the source has been considered by Quevedo 
and Mashhoon l26l. 
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For our study we refer to the solution given by Young and Coulter [19] , 
In the limit of weak fields and small quadrupole moments the nonvanish- 
ing Geroch-Hansen moments associated with this solution are the monopole 
A^o = M and the quadrupole moment A4 2 = Q, the latter being related to 
the parameter q by the equation 

w = y < ? > 

according to our conventions G = 1 = c. 

Due to the stationarity of this spacetime a suitable family of fiducial ob- 
servers is that of the so called static observers, with unit timelike four velocity 
e t - = e~^dt aligned with the timelike Killing vector dt- An orthonormal frame 
adapted to the static observers is given by 

e { - = e-^, e p = e^d p , e^ = -e^, e s = e^d z . (8) 



2.1 Circular orbits 

Circular orbits (p — const) onaz = const hypersurface have a four velocity 
of the form 

U = Tute + ue$), iu = (1 - v 2 )- 1 ' 2 , (9) 

with v Lie-constant along U, i.e. £tjv — 0. They form a 1-parameter family 
of orbits parametrized by v. In general these orbits are accelerated with a 
transverse acceleration a(U) = DU/dru — aPep + a z e£ such that 



p 



H=^e^{l + v 2 )i, z . (10) 



It is convenient to introduce a polar representation for the acceleration com- 
ponents, i.e. 

ap = KCOsx, af = Ksinx, (11) 

with 



n=ya 2 ~ + aj, tanx = a z /ap. (12) 

Therefore we have 

a(U) — n{co&x e p + sm X e z) = • (13) 

Starting a Frenet-Serret procedure with U = Eq one obtains a Frenet-Serret 
frame governed by the transport equations: 

DE Q DE x 

— — = kEi , — — = nE + t x E 2 , 

dm dru 

DE 2 DE 3 , , 

= - Tl E x + t 2 E z , — -i = -r 2 E 2 . 14 

dru dru 
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The curvature k is the magnitude ||a([/)|| of the acceleration a(U), while the 
first and second torsions T\ and Ti are the components of the Frenet-Serret 
angular velocity vector 

w (FS) = nE 3 + t 2 E x , ||w (F s)|| = [if + tI] 1/2 , (15) 
putting the spatial transport equations (|14l) in the form 
DE 

—j — - = w (FS) x E a + kE Si , (16) 

where x denotes ordinary vector product in the Euclidean three space or- 
thogonal to U. It is easy to show that in this case the Frenet-Serret frame 
vectors are 



Ei = cos xep + sin x^z , 

1 dU 

E 2 = M»e i + e $ ) = -^—, 
E 3 = - sin xep + cos %e. 



dEi 
dx 



while the Frenet-Serret torsions [27l[28] are given by 

1 dn k dx 



T\ = 



T2 = 



2jy dv ' dv 



(17) 



(18) 



Moreover, apart from those circular orbits located on the symmetry hyper- 
plane z — (where ip z = 0), no geodesies exist for whatever special choices 
of v. Vice versa, circular orbits on the symmetry plane z — correspond to 



implying also 



t 1 = 7 ^-7_(i_2^p) 



(19) 
(20) 



and T2 — 0. Note that in this case (symmetry hyperplane) we have fixed 
E\= ep and allowed k to vary its sign, as it is customary [28] . 

A direct evaluation shows that circular geodesies onz = correspond to 



(goo) 



± 



p^p 


1/2 




.1 - pi> P _ 




2=0 



(21) 



Using the solution dU)— (J4j) we find that 

~q-qo 



(goo) 



1/2 



(22) 
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where the quantities qo and q\ (which are functions of r*/Af) are the quadrupole 
critical values corresponding to pip p = and pip p = 1 respectively: 



9o 



3 7% 

4M 



_1L In 



M 2 



M 



M 



3^ 
2M 2 



1 



<7i = 9o 



1 - — 

M 



(23) 



Note that at a fixed value of p the quantity /3^ p is a linear function of q. 

The geodesic velocities (|22l) are plotted in Fig. [1] both as functions of the 
quadrupole parameter q for fixed radial distance (see Fig. (a)) and as func- 
tions of p/M for different values of q (see Fig. (b)). In the first case (Fig. (a)) 
we have shown how the quadrupole moment affects the causality condition: 
there exist a finite range of values of q wherein timelike circular geodesies 
are allowed (see the discussion below) . The difference from the Schwarzschild 
case is clear instead from Fig. (b): the behaviour of the velocities differs sig- 
nificantly at small distances from the source, whereas it is quite similar for 
large distances. 

For a later use it is useful to consider the following two limits of 



M/p — > for fixed q and q 
we have 



(geo) • 

for fixed M and p. In the first case {M/p — > 0) 



(gco) 



p 



1/2 



1 M 
1+r I — 
P 



In the second case (q — ► 0) we have instead 

"(to) = ±v k ± qW{v K ) 

where 

1 



4g-5 
10 



0(q 2 ) , 



7/2 



(24) 



(25) 



8z4 



* [3(1 + v\){\ + 2v\) ln(l + 2v%) - 2i%(v% + + 3)] 



W{v K ) 

(26) 

is a function of the linear velocity vk of circular geodesies in the Schwarzschild 
spacctime given by 



vk 



M 



It is easy to show that in the limit M/p — > one has 

1/2 



(27) 



v K — 



1 / M 
1+r I — 

P 



5/2 



1 



7/2 



01^ 



7/2 



(28) 



Let us discuss now the causal properties of geodesies. The requirement 
that the argument of the square root of Eq. (|2"Tj) be positive gives 



P^p 
1 - P^p 



> 



< pn\) p < 1 



qi <q<q Q 



(29) 
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(a) (b) 

Fig. 1 The geodesic linear velocities ^ eo N are plotted in Fig. (a) as functions of the 
quadrupole parameter q for fixed radial distance p/ M = 4 from the source. Time- 
like circular geodesies exist for q2 < q < qo (see Eq. (I31|l 'l. with go ~ 82.84 and qi ~ 
—87.94 for this choice of parameters. The behaviour of v f geo \ as a function of p/M is 

shown in Fig. (b) for different values of q = [-80, -30, -5, 0, 2.5, 5, 10, 20, 50]. The 
case q — (Schwarzschild) corresponds to the dashed curve without a relative max- 
imum. The shape of the curves with negative q is very similar to the Schwarzschild 
one, the lightlike condition = 1 being reached at greater values of the radial 

distance for decreasing values of q. The curves with positive q instead all present 
a relative maximum and are ordered from left to right for increasing values of q. 
Shaded regions are forbidden. 



where the quantities qo and q± are given by Eq. (|23[) . The further requirement 
that |f^g eo \| < 1 ensures that they remain timelike. This condition implies 

<? > 7^(9o + qi) = <72 , (30) 

so that the condition for the existence of timelike circular geodesies turns out 
to be 

92 < q < qo ■ (31) 

Fig. [2] shows the behaviours of the critical quadrupole parameters qo, qi and 
qi as functions of the radial distance p/M. 

In terms of the quantities qo and q\ the expression (|19p for the acceleration 
k becomes 
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100 



qO 



5 



6 



7 



p/M 



q2 



-100 



Fig. 2 The behaviours of the critical quadrupole parameters go, qi and qi are 
shown as functions of p/M. 

Its behaviour as a function of v is shown in Fig. [3] for different values of the 
quadrupole parameter. For increasing values of q the corresponding values 
of ^ eo ) decrease up to a critical value (corresponding to = 0) beyond 

which they do not exist anymore. 

A similar discussion can be done for the first torsion, i.e. there exists 
a critical value of the quadrupole moment q such that n = for any v. 
This happens when pij) p = 1/2. It is quite surprising that the corresponding 
critical value of q coincides with the limiting value qi for timelike geodesies 
as defined in Eq. (f30| . In fact, by introducing the quantities qo, qi and q<± 
Eq. (|20|) can be cast in the form 



We thus recognize how the presence of the quadrupole moment q changes 
completely and enriches the situation with respect to the Schwarzschild case, 
examined long ago by many authors [27l[28] . This is well elucidated in Fig. [4] 
where the behaviour of t± as a function of v is discussed for different values 
of the quadrupole parameter. 

3 Motion of quadrupolar particles in the Schwarzschild 
background 

We consider now the complementary case of an extended body with structure 
up to the quadrupole in the field of a Schwarzschild background, following the 




(33) 
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-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.! 



(a) 



(b) 



Fig. 3 The magnitude of the acceleration k for circular orbits at z = is plotted as 
a function of v for p/M = 4 and different values of the quadrupole parameter: (a) 
q = [0, 20, 40, 60, 80, 100] and (b) q = [100, 200, 300, 400, 500, 600]. The (symmetric) 
values of v associated with k = correspond to geodesies, i.e. ^ eo N (see Fig. (a)). 
The various curves are ordered so that for increasing values of q the value of l^^,) | 
decreases until it approaches the critical value go ~ 82.84 beyond which circular 
geodesies do not exist anymore (see Fig. (b)). It is worth noting that the situation 
is very similar for selected negative values of q. 



description due to Dixon p~2lfT3irT4i ri 5.16]. In the quadrupole approximation 
Dixon's equations ar^3 



DPf 



^obU" 'S ap - \j a(SlS R a p-,s^ = F (spin)M + F^ naA > ,(34) 
dry 2 6 



dru 



= 2P^U u] + - J^^ifl a = 2P^U v] + D^^^ , 



(35) 



where P^ = mU£ (with U p ■ U p = —1) is the total four-momentum of the 
particle, and is a (antisymmetric) spin tensor; U is the timelike unit 
tangent vector of the "center of mass line" Cjj used to make the multipole 
reduction, parametrized by the proper time tjj- The tensor J a P"< s is the 
quadrupole moment of the stress-energy tensor of the body, and has the same 

1 Note that the torque t erm p( quad ) ^ " in the second set of equations was mis- 
printed in previous works 29, 30,31 32]. The results of all those papers are but 
correct. 
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(a) (b) 

Fig. 4 The first torsion n for circular orbits at z = is plotted as 
a function of v for p/M = 4 and different values of the quadrupole 
parameter: (a) q = [-600,-500,-400,-300,-200,-100] and (b) q = 
[-80,-40,0,100,200,300,400,500,600]. The behaviour changes depending on 
whether q is less or greater than the critical value §2 ~ —87.94. The curves in 
both Figs, (a) and (b) are ordered in such way that for increasing \q\ they shrink 
on the vertical axis. 



algebraic symmetries as the Riemann tensor, i.e. 20 independent components. 
Note that there are several equivalent expressions used in the literature for 
the torque £)(i uad )A« / j which can be summarized by 

n (quad) M „ _ 4 r TV] a/37 n(l uacl )/"' — _ f>v] /o fi \ 

•^dixon — £ n a01 J ! LJ \ 1 -t _ 3 lafi ' ^ > 

due to Dixon ( 16], p. 65) and Ehlers and Rudolph ([33], p. 209). 

There are no evolution equations for the quadrupole as well as higher 
multipoles as a consequence of the Dixon's construction, so the structure 
only depends on the considered body. Therefore the system of equations is 
not self-consistent, and one must assume that all unspecified quantities are 
known as intrinsic properties of the matter under consideration. 

Moreover, in order the model to be mathematically correct certain ad- 
ditional conditions should be imposed [T2]. As it is standard one may limit 
considerations to Dixon's model under the further simplifying assumption 
34,33] that the only contribution to the complete quadrupole moment J a Pi s 
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stems from the mass quadrupole moment Q a @ so that 

J a ^ s = SU^Q^U^ , Q af) U pP = . (37) 

We are interested here in bodies with vanishing spinning structure, so 
that the model equations undergo big simplifications and reduce to 



Dm i 



in- 



Atu = 2 U P Q07U P Ra W" ' (38) 
mUplTl = U«Q^U]?I?\ a p - U«U;Q^R v \ a p , (39) 

where we have also assumed for simplicity the mass of the body as a constant. 

Let us consider the case of the Schwarzschild spacetime, with the metric 
given by Eq. |T]) with 7 = 70 and if> = ipo , given by Eq. . 

Let us assume that U is tangent to a (timelike) spatially circular orbit 

U = 7u [ ei + ue $ ], <yt, - (1 - v 2 )- 1 ' 2 , (40) 

with v constant along U. We limit our analysis to the equatorial plane [z = 0) 
of the Schwarzschild solution. The linear velocity corresponding to circular 
geodesies is given by Eq. ([27]) , whereas the angular velocity is given by 



r = — = v * 

Let also P = m U p be such that 

U p = 7 p h + v v e^\ , 7p = (1 " ^y 1 ' 2 . ( 42 ) 

i.e. let us assume that U p also is tangent to a circular orbit and set up an 
orthonormal frame adapted to U p given by 

e = U P , ei = e p , e 2 = e~ z , e 3 = j p \v p e~ t + e^] ; (43) 

hereafter all frame components of the various fields are meant to be referred 
to such a frame. Note that the assumption of having both U and U p leads to 
great simplifications to Dixon's equations. 
From Eq. ([37)> 9 we have 



/on 



= Q01 = Q02 = O03 = . (44) 



We also assume that that all the surviving components of the mass quadrupole 
moment are all constant along the path. The latter assumption corresponds 
to the definition of "quasirigid motion" (or "quasirigid bodies" ) due to Ehlers 
and Rudolph |S3] . Clearly in a more realistic situation the latter hypothesis 
should be released. 

Consider first the constraint equations (f3"9")) . They imply that 



Q12 = Q13 = Q23 — , 



(45) 
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and after introducing the following "structure functions" of the extended 
body 

Qn = Q22 - f, Q33 = Q22 - f , (46) 

they also give 

= 7£/ (z,-z, p )-3 7p ^/ . (47) 
m 

Note that the quantities / and /' are necessarily small, in order to avoid 
backreaction effects. 

Consider then the equations of motion (|38[) . They imply that 

= ju7 P (w v - v\) - l^[f'+ (1 - 3 7p 2 )/] ■ (48) 

Solving Eqs. (fiT|) and (I48[) for v and v v in terms of/, /' completely determines 
the motion. 

The quadrupolar structure of the body turns out to be represented by 
three quantities: /, /' and €£22- However, classically, the quadrupole moment 
tensor of a mass distribution is tracefree. Assuming for simplicity the same 
property to hold also for the relativistic quadrupole moment tensor implies 

= Qn + Q22 + Q33 = 3Q 22 -/-/', (49) 

so that the components Q a b in this case are completely determined by / and 
/' only 

Qxx = -\f + \f, Qn = \tf + f), Q™ = \f-\f ■ (50) 

If the body is axially symmetric about the z-axis, then f' = f and the frame 
components of Q reduce to 

Q o6 = diag[-//3,2//3,-//3] . (51) 

It is worth noting that the above assumptions for the extended body's struc- 
ture (i.e. constant frame components for the quadrupole tensor represented 
by a single structure function /) here used in order to make our analysis as 
simple as possible can be easily relaxed in favour of a more realistic descrip- 
tion. 

Eqs. (|47]) and j48]) then become 

f 

= 717(1/ - u p ) - ZlpV P C,K— > ( 52 ) 

- ful P {yu p - v\) - I (2 - 3 7p 2 )Ck^ • (53) 

The above relations define the kincmatical conditions allowing circular mo- 
tion of the extended body taking into account its quadrupolar structures. 

Due to the smallness of quadrupolar quantity / the motion of the ex- 
tended body will be nearly geodesic with correction which we will retain up 
to first order in /. Eqs. and (fB"3"]) thus imply 

vf> ~v ± t*vk&- , (54) 



VK 



3Ci 



13 



where the signs ± correspond to co/counter rotating orbits. The correspond- 
ing angular velocity £± = /vk)v± and its reciprocal are 



1 



4*4 



1 



1 

±5 



1 



3&/ 



4-1 



(55) 



4 Quadrupolar particles interacting with gravitational monopoles 

The linear velocity (f5"4"]h of the center of mass of the extended body described 
using Dixon's model, namely 



v ± ~± 



3&/ 



4z/ft- TO 



(56) 



should now be compared with the geodesic linear velocities -J co ) given by 

Eq. (|25|) . in the limits in which comparison is really allowed (i.e. in the limit 
of validity of Dixon's model): small mass and quadrupole moment of the 
body if compared with the mass of the central object. Note that a similar 
treatment was done in the case of a spinning particle orbiting a Schwarzschild 
black hole compared with a geodesic in the Kerr spacetime [55] , 
The result (to first order in both q and /) is 

-T* L £; = QW{vk). (57) 
Avk to 

In the limit of large distances from the central source (p 3> M) Eq. ([57)) gives 
the nice result 

L = l M 2 q = % . (58) 
to 15 M v ' 

This correspondence validates once more Dixon's model. 



5 Concluding remarks 

We have considered the static gravitational field of a quasi-spherical source 
belonging to the Weyl class of solutions of Einstein's equations. In this field 
we have discussed the geometric properties of (accelerated in general) time- 
like circular orbits, by analyzing the associated Frenet-Serret curvature and 
torsions. Among these orbits we have studied in detail the geodesies on the 
equatorial plane and found a number of interesting properties. For instance 
we have found that there exists a critical value for the quadrupole moment 
beyond which the are no more timelike circular geodesies. This fact was unex- 
pected, on the basis of previous works concerning static Schwarzschild black 
hole spacetime. 

However, the most important contribution of this paper concerns com- 
parison between geodesic motion in the field of a quasi-spherical source with 
the motion of an extended body (quadrupolar particle, described by using 
Dixon's model) in the background of a Schwarzschild black hole. We have 
shown how the quadrupole moment of a gravitational source as described 
by the Geroch-Hansen approach coincides with the quadrupole moment of a 
Dixon's extended body. 
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